6034 





NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


BEPORT No. 865 

METHOD FOR CALCDLATmC WmG OT ARACTERISTICS 
BY LIFTING-LINE THEORY USING NONLINEAR 
SECTION LIFT DATA 


By JAMES C. SlYELLS and ROBERT H. NEELY 



1947 


For Mie br UMSnperlntendeiit of Documento, U.S. GoTemment Printing Office, Washington 25. D. C. ......... Price M cents 
















w 

Q 

m 

I 

M 

8 

a 

b 

c 

A 

V- 

3 

L 

D 

D, 

Dy 

D, 

a 


AERONAUTIC SYMBOLS 


1. FUNDAMENTAL AND DERIVED UNITS 



Syoobol 

Metric 

YlngKah 

Unit 

Abbre-vlar 

tion 

- ; Unit - ^ 

Abbrevia- 

tion 


1 

t 

F 

niAf«n> 

/ m 

S' ■ 

foot (or mile)-. 

ft (or mi) 
eec (or hr) 
lb 


fWAnrifl 

second (or. hour) — 

Force 

weight of 1 kilogram-...— 


weight of 1 pound- 

\ ' 

P 

V 

Vinrw^y*W**r (nrifitrift) . . 


horsepower - 

hp 

mph 

fps 

Bpeed-- — -i- 

ridlometers per hour 

\meters per seCond 

- ]^h 
. mps. 

miles'per hour. — 

feet i>er second — 


2. GENERAL SYMBOLS 


Wdght=*nff _ ^ 

Standard acceleration of gravity=9.80665 m/s* 
or 32.1740 ft/sec* 

Mass=~ , ^ 

Moment of inertia—wi*. Qbidicate axis of 
radius of gyration k by proper subksript.) 
Coefficient of viscosity 


p ■ Kinematic viscosity 
p . Density (mass per unit volume) 

Sfemdard density of dry air, 0.12497 at 16® O 

and 760 mm; or 0.002378 lb-ft“* sec* 

Specific weight of “stand^d” air, 1.2255 kg/in* or 
0.07651 Ib/ou ft ,, : ■. 


8. AERODYNAMIC SYMBOLS L 7 


Area 

Area of wing 


Spa,n ' 

Chord 

Asp^t ratio, ^ 

True air speed 
Dynamic pressure, 


Lift, absolute coefficient Cl=-^ 


D 


Drag, absolute coefficient Cd=-^ 
Profile drag, absolute coefficient 


Do 


D, 


Induced drag, absolute coefficient CD^=■^ 




Parasite drag, absolute coefficient 
Cross-wind force, absolute coeffici^t 


C 


iic 

it 

Q 

Q 

B 


a 

e 

ao 

Oti 

Ola 



Angi ft of setting of wings (relative to thrust line) 
Angle of stabilizOT setting (relative to thrust 
Ime) 1 

Resultant moment ^ v , 

Resultant angular velocity .. 

Reynolds number, p^ where lis a linear dimen- 
sion (e.g., for an airfoil of 1 .0 ft chord, 100 mph.7 
. standard pressure at. 15° C> the corresponding 
Reynolds number is 935,400; or for an airfoil 
of 1.0 m chord, 100 mps, the corresponding 
Reynolds nmnber is 6,866,000) 

Angle of attack - 

Angle of downwash 

Angie of attack, infinite aspect ratio . - 

Angle of attack, induced 

Angle of attack, absolute (measure from zero- 
lut position) 

Plight-path angle 
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METHOD FOR CALCULATING WING CHARACTERISTICS BY LIFTING-LINE THEORY USING 

NONLINEAR SECTION LIFT DATA 

By James C. Sivells and Robert H. Neely 


SUMMARY 

A method is ■presented jor calculating wing characteristics by 
lifting-line theory using nonlinear section lift data. Material 
from various sources is combined with some original work into 
the single complete method described. Multhopp’s systems of 
multipliers are employed to obtain the induced angle of attack 
directly from the spanwise lift distribution. Equations are 
developed for obtaining these multipliers for any even number 
of spanwise stations, and values are tabulated for 10 stations 
along the semispan for asymmetrical, symmetrical, and anti- 
symmetrical lift distributions. In order to minimize the com- 
puting time and to illustrate the procedures involved, simplified 
computing forms containing detailed examples are given for 
symmetrical lift distributions. Similar forms for asymmetrical 
and antisymmetrical lift distributions, although not shown, can 
he readily constructed in the same manner as those given. The 
adaptation of the method for use with linear section lift data 
is also illustrated. This adaptation has been found to require 
less computing time than most existing methods. 

The wing characteristics calculated from general nonlinear 
section lift data have been found to agree much closer with ex- 
perimental data in the region of maximum lift coefficient than 
those calculated on the assumption of linear section lift curves. 
The calculations are subject to the limitations of lifting-line 
theory and should not be expected to give accurate results for 
wings of low aspect ratio and large amounts of sweep. 

INTRODUCTION 

The lifting-line theory is the best known and most readily 
applied theory for obtaining the spanwise lift distribution of 
a wing and the subsequent determination of the aerodynamic 
characteristics of the wing from two-dimensional airfoil data. 
The characteristics so determined are in fairly close agree- 
ment with experimental results for wings with small amounts 
of sweep and with moderate to high values of aspect ratio; 
for this reason, this theory has served as the basis for a large 
part of present aeronautical loiowledge. 

The hypothesis upon which the theory is based is that a 
lifting wing can be replaced by a lifting line and that the 
incremental vortices shed along the span trail behind the 
wing in straight lines in the direction of the free-stream 
velocity. The strength of these trailing vortices -is propor- 
tional to the rate of change of the lift along the span. The 
trailing vortices induce a velocity normal to the direction of 
the free-stream velocity and to the lifting line. The effec- 
tive angle of attack of each section of the wing is therefore 


different from the geometric angle of attack by the amount 
of the angle (called the induced angle of attack) whose 
tangent is the ratio of the value of the induced velocity at ; 
the lifting line to the value of the free-stream velocity. The 
effective angle of attack is thus related to the lift distribution 
through the induced angle of attack. In addition, the ef- 
fective angle of attack is related to the section lift coefficient 
according to two-dimensional data for the airfoil sections in- 
corporated in the wing. Both relationships must be simul- 
taneously satisfied in the calculation of the lift distribution 
of the wing. 

If the section lift curves are linear, these relationships may 
be expressed by a single equation which can be solved analyt- 
ically. In general, however, the section lift curves are not 
linear, particularly at high angles of attack, and analytical 
solutions are not feasible. The method of calculating the 
spanwise lift distribution using nonlinear section lift data 
thus becomes one of making successive approximations of 
the lift distribution unitl one is found that simultaneously 
satisfies the aforementioned relationships. 

Such a method has been used by Wieselsberger (reference 1) 
for the region of maximum lift coefficient and by Boshar 
(reference 2) for high-subsonic speeds. Both of these writers 
used Tani’s system of multipliers for obtaining the induced 
angle of attack at five stations along the semispan of the 
wing (reference 3). Tani, however, considered only the case 
of wings with symmetrical lift distributions. Multhopp 
(reference 4), using a somewhat different mathematical treat- 
ment from that which Tani used, derived systems of multi- 
pliers for symmetrical, antisymmetrical, and asymmetrical 
lift distributions for 4, 8, and 16 stations along the semispan. 
Multhopp’s derivation, in slightly different form and nomen- 
clature, is presented herein and tables are given for the 
multipliers for 10 stations along the semispan (the usual 
number of stations considered in many reports in the United 
States). 

For synimetrical distnbutiohs of wing chord and angle of 
attack, the multipliers for symmetrical lift distributions may 
be used with nonlinear or linear section lift curves. For 
asymmetrical distributions of angle of attack, the multi- 
pliers for asymmetrical lift distributions must be used if 
nonlinear section lift curves are used. If an asymmetrical 
distribution of angle of attack can be broken up into a 
symmetrical and an antisymmetrical distribution, the anti- 
symmetrical part may be treated separately if the section 
lift curves can be assumed to be linear. 
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The purpose of the present paper is to combine the con- 
tributions of Multhopp and several other writers, together 
with some original work, into a single complete method of 
calculating the lift distributions and force and moment 
characteristics of wings, using nonlinear section lift data. 
Simplified computing forms are given for the calculation of 
symmetrical lift distributions and their use is illustrated by 
a detailed example. The adaptation of the method for use 
with linear section lift data is also illustrated. No forms 
are given for asymmetrical or antisymmetrical lift distri- 
butions inasmuch as such forms would be very similar to 
those given. 
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SYMBOLS 


wing area 
wing span 

chord at any section 
root chord 
tip chord 

mean geometric chord iSIh) 

( 2 

■g I cHy 


aspect ratio {h^jS) 
coordinate parallel to root chord 
coordinate perpendicular to plane of symmetry 
coordinate perpendicular to root chord and parallel 
to plane of symmetry 

free-stream dynamic pressure pV^ 


Reynolds number {pVcIfi or pVc' jy) 
mass density 
free-stream velocity 
coefficient of viscosity 
wing lift coefficient (LfqS) 
section lift coefficient (l/qc) 
wing lift 
section lift 

wing drag coefficient (D/qS) . 
wing profile-drag coefficient 
wing induced-drag coefficient 
section profile-drag coefficient 
section induced-drag coefficient 
wing drag 

wing pitching-moment coefficient (M/qSc') 
section pitching-moment coefficient about section 
quarter-chord point 
wing pitching moment 
wing rolling-moment coefficient {L' /qSb) 
wing rolling moment 

wing induced-yawing-moment coefficient 
wing profile-yawing-moment coefficient 
angle of attack of any section along the span 
referred to its chord line 

angle of attack of root section referred to its chord 
line 

angle of attack of root section referred to its zero 
lift line 


at 

ao 



®S{Z,=0) 

e 


/ 

Ct 

a 

ao 


cos 8 



ymk 

Vm 

Vms 

^ma 

E 


section induced angle of attack 
effective angle of attack of any section 
section angle of attack for two-dimensional airfoils 
angle of zero lift of any section 
angle of zero lift of root section 
wing angle of attack for zero lift 
geometric angle of twist of any section along the 
span (negative if washout) 
aerodynamic angle of twist of any section along the 
span (negative if washout) 
geometric angle of twist of tip section 
aerodynamic angle of twist of tip section 
wing lift-curve slope, per degree 
section lift-curve slope, per degree 
/ Two-dimensional lift-curve slope \ 

V Edge-velocity factor / 
coordinate (2ylb) 
coefficients in trigonometric series 
multiplier for induced angle of attack (asymmetrical 
distributions) 

multiplier for induced angle of attack (symmetrical 
distributions) 

multiplier for induced angle of attack (antisym- 
metrical distributions) 

multiplier for lift, drag, and pitching-moment 
coefficients (asymmetrical distributions) 
multiplier for lift, drag, and pitching-moment 
coefficients (symmetrical distributions) 
multiplier for rolling- and yawing-moment coeffi- 
cients (asymmetrical distributions) 
multiplier for rolling-moment coefficient (anti- 
symmetrical distributions) 

edge-velocity factor J 


Subscripts: 

max ridaximum value 

al value for additional lift ((7 l=1) 

b value for basic lift (C'^=0) 

(ttoj) value for constant value of 

(e/) value for given value of «/ 


THEORETICAL DEVELOPMENT OF METHOD 


LIFT DISTRIBUTION 


The methods of Tani (reference 3) and Multhopp (refer- 
ence 4) for determining the induced angle of attack are 
fundamentally the same, differing only in the mathematical 
treatment. The method presented herein is essentially the 
same as that given by Multhopp. In the following deriva- 
tion the spanwise lift distribution is expressed as the trig- 
onometric series 


^-^=J^AnSmnd ( 1 ) 

2v 

as in reference 5, where 8 is defined by the relation cos 6— 

It may be noted that each coefficient An, as used herein, is 
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equal to four times the corresponding coefficient in reference 5. 
The induced angle of attack (in degrees) at a point yx on the 
lifting line is 

. TT 87rJ-fe/2 dy 
Vi—y 


■dy 


( 2 ) 


This integral (in different nomenclature) was given by 
Prandtl in reference 6. If equation (1) is substituted into 
equation (2) and the variable is changed from y to B, the 
induced angle of attack at the general point B becomes, 
according to reference 5, 


180 ^ ^ ^ 


(3) 


The problem of obtaining the induced angle of attack is thus 
reduced to one of determining the coefficients of the trigo- 
nometric series. 

The lift distribution (equation (1)) may be approximated 
by a finite trigonometric series of r — 1 terms where, for sub- 
sequent usage, r is assumed to be even. The values of Cicjh 

at the equally spaced points in the range 0<^<7 t are 

expressed as 


\ " / m 71=1 


Sin n 


mw 


(4) 


where ?n=l, 2, 3, . . . r — 1. Conversely, if the values of 
Cjc/6 are known at each point, the coefficients An of the finite 
series may be found by harmonic analysis as 


sm-a 

' 7n = l V ^ /m 


rr^ 

r 


(5) 


If equation (5) is substituted in equation (3), a double 
summation is obtained for the induced angle of attack as 

180 /c;c\ . rmr~\ 

4tt sin 0 5^1 V 0 Jm r J 


47rr sin0 


~JA 


If the induced angle of attack is to be determined at the 
same points 0 at which the load distribution is known, that 
kir 

is, at the points 0=—; then 


180 


** . . kir 

4rr sin — 

= Pm. 

77i = l\ 0 /m 


s(^) s»r 


(k — m)Tr (Ar-|-m)x 

cosa^^ cosa^= — 


] 


( 6 ) 


where 


180 


47rr sin 


kv 


r-l r 

' 

7» = 1 L 


(k — Tn,)ir 
cos n~ cos n 




It can be shown that, if cos 

rcos(r— 1)^ — (r — l)cosr^ — 1 


„=i 2(1— cos 0) 


^ n cos ruf> 

71 — I 

If ^=0, a numerical series is obtained 

n=l ^ 

By use of these relationships in equation (7) it is found that, 
when is odd, 


Pmk~ 


180 


. . k-K 

4-Trr sm — 
r 


1 


, (k+m)TT , {k—m)tc 

1 — cos^^ — i — — 1— cos^^ ^ 


when k=m, 


P. 


180r 


. kir 
87t sm — 


(8a) 


(8b) 


(8c) 


and when ^±m is even and kr^yn, 

Pm. = Q 

For a symmetrical lift distribution 

V b A'V b 

and 


so that the summation for needs to be made only from 1 


to r/2 


■»<,=f:(-f) 

77i=l \ 0 / 71 

where, when is odd, 




(9) 


^mk — Pmk~^ Pr—m, k 

_ 180 

2irr sm — 
r 


(for 


r r 


. (Ar+m)TT 


sin 


sin 


(fc— m)x 


(lOa) 


^mS; Pmk ^for 7fl 2 ^ 


180 


/ 2fcir , A 
Trr (^cos — -rl J 


(10b) 


when k=m, 


^mk — Pmk 


180r 


„ • kir 

8tt sm — 
r 


and when is even and kj^vHy 

— 0 


(10c) 


(lOd) 
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For an antisymmetrical lift distribution 

\ b Jm \ ^ / T—m 

and 

In this case the summation for ai^ needs to be made only 
from 1 to 1 since =0; then 


^ \ ~h~ ) 

m — l\ 0 


ymk 


( 11 ) 


where, when k±?n is odd, 

T mb ^mk ^T—m, k 


2irt' 


1 


1 




sim 


sm^ 


(12a) 


when k—m, 


ymk — ^mk 

180r 


. kir 
8tt sin — 
r 


(12b) 


and when is even and k^m, 

7m/t = 0 


(12c) 


Multipliers can thus be calculated so that the induced 
angle may be readily obtained by multiplying the known 
values of Cic/b by the appropriate multipliers and adding the 
resulting products. The multipliers are independent of the 
aspect ratio and taper ratio of the wing. Tables I and II 
present values of ^mk, and \mk and jmk, respectively, for 


r=20. 


Similar tables for 


4ir . . 4 tt . . 

and ai e given m 


TABLE I.— INDUCED-ANGLE-OF-ATTACK MULTIPLIERS FOR ASYMMETRICAL LIFT DISTRIBUTIONS ■ 


[“’‘"JiCtO./'"*] 


\ 

2k 

b 













\ 

\ 

\ 

-0. 9877 

-0.9511 

-0. 8910 

-0. 8090 

-0. 7071 

-0. 5878 

-0, 4540 

-0. 3090 

-0. 1564 

0 



21/ 

b 

* 

\ 

19 

18 

17 

16 

15 

14 

13 

12 

11 

10 



-0. 9877 

19 

915 051 

—106. 985 

0 

-7. 019 

0 

-1.401 

0 

-0. 486 

0 

-0. 230 

1 

0. 9877 

-.9511 

18 

-329. 859 

463. 533 

-122. 749 

0 

-7. 438 

0 

-1.792 

0 

-.701 

0 

2 

.9511 

-.8910 

17 

0 

-180. 336 

315. 512 

-96. 737 

0 

—7. 073 

0 

-1.920 

0 

-.819 

3 

.8910 

-.8090 

16 

-26. 374 

0 

-125.246 

243. 694 

-81. 067 

0 

-6. 680 

0 

-1.977 

0 

4 

.8090 

-. 7071 

15 

0 

-17.020 

0 

-97. 524 

202. 571 

-71.139 

0 

-6. 391 

0 

-2. 026 

5 

.7071 

-.5878 

14 

-7.246 

0 

-12.604 

0 

-81. .392 

177.054 

-64. 735 

0 

-6. 228 

0 

6 

. 5878 

-.4540 

13 

0 

-5. 166 

0 

-10.126 

0 

-71.296 

160. 761 

-60. 725 

0 

-6. 192 

7 

. 4540 

-.3090 

12 

-2.958 

0 

-4. 022 

0 

-8. 596 

0 

-64. 817 

150. 611 

-58. 514 

0 

8 

. 3090 

-. 1564 

11 

0 

-2. 241 

0 

-3. 322 

0 

-7. 004 

0 

-60. 768 

145. 025 

-57. 812 

9 

. 1.564 

0 

10 

-1. 468 

0 

-1.804 

0 

-2. 865 

0 

-6. 950 

0 

-58. 533 

143. 239 

10 

0 

.1564 

9 

0 

-1.153 


-1.518 

0 

-2. 554 

0 

-6. 530 

0 

-57. 812 

11 

-. 1.564 




0 

-.946 

0 

-1.319 

0 

-2.340 

0 

-6. 288 


12 

-.3090 




-.646 



0 

-1.176 

0 

-2. 192 

0 


13 

-. 4540 

.5878 

6 

-.467 

0 


0 

-.691 

0 

-1.068 



0 

14 

-.5878 

,7071 

6 

0 

-..368 

0 

-.441 


-.604 

0 



-2. 026 

15 

-. 7071 


4 

-.261 


-.291 


-.366 

0 

-.528 

0 


0 

16 

-.8090 

.8910 

.3 

0 

-. 192 

0 

-.225 

0 

-.297 

0 

-.452 

0 

-.819 

17 

-.8910 

. Qr'in 

2 

-. 118 

0 

-. 130 

0 

-. 161 

0 

-.224 


-.361 

0 

18 

-.9.511 

.9877 

1 

0 

-.060 

0 

-.069 

0 

-.090 

0 

-. 133 

0 

-.2-30 

19 

-.9877 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

^ \ 

2.1/ 

b 



.9877 

.9511 

.8910 

.8090 

.7071 

. 5878 

.4540 

.3090 

.1564 

0 

• 2?/ 
6 

\ 

\i 


I Values of fc at top to be used with values of m at left side; values of k at bottom to be used with values of m at right side. 
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TABLE IL— INDUCED-ANGLE-OF-ATTACK MULTIPLIERS FOR SYMMETRICAL LIFT DISTRIBUTIONS AND FOR 

ANTISYMMETRICAL LIFT DISTRIBUTIONS 



0 

0. 1564 

0.3090 

0. 4540 

0. 5878 

0. 7071 

0. 8090 

0. 8910 

0. 9511 

0. 9877 

b 

TO\ * 

10 

9 

8 

7 

6 

5 

4 

3 

2 

1 

10 /C(C\ 

Multipliers Xmt S ( ) Xmt 

0 

10 

143. 239 

-58.533 

0 

-6.950 

0 

-2. 865 

0 

-1.804 

0 

-1. 468 

.1564 

9 

-115. 624 

145. 025 

-07. 298 

. 0 

-10.158 

0 

-4. 840 

0 

-3.394 

0 

.3090 

8 

0 

-64. 802 

150.611 

-67. 157 

0 

-9. 916 

0 

-4.968 

0 

-3.768 

.4540 

7 

-12.384 

0 

-62. 917 

160. 761 

-72.472 

0 

-10. 920 

0 

-5.812 

0 

.5878 

6 

0 

-8. 320 

0 

-65. 803 

177. 054 

-82. 083 

0 

-13. 134 

0 

-7.713 

.7071 

5 

-4.051 

0 

-7.372 

0 

-71. 743 

202. 571 

-97.965 

0 

-17.388 

0 

.8090 

4 

0 

-2. 880 

0 

-7.208 

0 

-81. 434 

243. 694 

-125. 537 

0 

-26. 035 

.8910 

3 

-1.638 

0 

-2.371 

0 

-7. .370 

0 

-96.902 

315.512 1 -180.528 

0 

.9.511 

2 

0 

-1. 062 

0 

-2.016 

0 

-7.599 

0 

-122.880 1 463.533 

-329.976 

.9877 

1 

-. 4.59 

0 

-.020 

0 

-1.491 

0 

-7. 089 

0 1 -167.045 

915. 651 

9 / \ 

Multipliers 7m* oti^= 

0. 1564 
. 3090 

9 

1 145. 025 

-54. 237 

0 

-5. 049 

0 

-1.804 

0 1 -1.087 

0 

8 

1 -.52. 226 

150. 611 

-62. 477 

0 

-7. 277 

0 

-3. 076 1 0 

-2. 147 

.4540 

7 


0 

-58.533 

160. 761 

-70. 120 

0 

-9. 326 

0 

-4.519 

0 

.5878 

r» 


-4. 136 

0 

-63. ()68 

177. 0.54 

-80. 701 

0 

-12.074 

0 

-6. 779 

.7071 

5 


0 

-5.410 

0 

-70. 535 

202. 571 

-97. 084 

0 

-16. 651 

0 

. 8090 

4 

1 

-1.074 

0 

-6. 152 

0 

-80. 701 

243. 694 

-124.955 

0 

-26. 113 

. 8910 

3 

1 

0 

-1.468 

0 

-6. 775 

0 

-96.512 

315. 512 

- 180. 145 

0 

.9511 

2 


-.340 

0 

-1.567 

0 

-7. 277 

0 

-122.619 

463. 533 

-329. 741 

. 9877 1 1 j 

0 

-.3.53 

0 

-1.311 

0 

-6.9.50 

{} 

-166. 926 

915. 651 


roferencc's 7 and 8, respectively, but no derivation is given 

therein. Tables for K,k, and Tmt are given in 

reference 4 for values of r—S, 16, and 32. An inspection of 
tables I and II shows that positive values occur only on the 
diagonal from upper left to lower right and that almost half 
of the values are equal to zero. The multipliers and X,„* 
may be used with either nonlinear or linear section lift data, 
whereas the multipliers for may be used only with linear 
section lift data. 

The method of determining the lift distribution becomes 
one of successive approximations. For a given geometric 
angle of attack, a distribution of Cj is assumed from which 
the load distribution Cic/b is obtained. The induced angle 
of attack is then determined by equation (6), (9), or (11) 
through the use of the appropriate multipliers and sub- 
tracted from the geometric angle of attack to give the 
effective angle of attack at each spanwise station. From 


section data for the appropriate airfoil section and local 
Reynolds number, values of Ci are read which correspond 
to the effective angle of attack of each section. If these 
values of C; do not agree with those originally assumed, a 
second assumption is made for Ci and the process is repeated. 
Further assumptions are made until the assumed values of 
Cl are in agreement with those obtained from the section 
data. 

WING CHARACTERISTICS 

Once the lift distribution of a wing has been determined, 
the main part of the problem of calculating the wing char- 
acteristics is completed. The induced-drag and induced-, 
yawing-moment coefficients are entirely dependent upon 
the lift distribution and it is assumed that the section profile- 
drag and pitching-moment coefficients are the same func- 
tions of the lift coefficient at each section of the wing as 
those determined in two-dimensional tests. 
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The calculation of each of the wing coefficients involves 
a spanwise integration of the distribution of a particular 

function J integration can be performed numeri- 

cally through the use of additional sets of multipliers which 
are found in the following manner. 


If 

/^^^=/(cos 6)='2fAn sin nO 

then 

J* f sin nd) sin 6 dd 


where 


7T . 2m7T 
^77l~~0 SlITl 
TO 8r r 


If the distribution is antisymmetrical,/^^^ = — 


where 

(^ma — 2(Tj)j 

Values of Vm, Vms, and are given in table III for r=20. 
TABLE III.— WING-COEFFICIENT MULTIPLIERS 


— A 
-2^1 


Since the values of are determined at the points 

Ai can be found by harmonic analysis as in equation (5) 


A 


_2 /2y\ 


sin 


m-ir 


Therefore 




2j/ 

b 

m 

Vm 



O’m a 

-0. 9877 

19 

0. 01229 


-0. 00607 


-.9511 

18 

. 02427 


-.01154 


-.8910 

17 

.03066 


-.01589 


-.8090 

16 

. 04616 


-. 01867 


-, 7071 

15 

. 05554 


-.01964 


-.5878 

14 

.06354 


-.01867 


-. 4540 

13 

.06998 


-.01589 


-.3090 

12 

. 07470 


-.01154 


-. 1564 

11 

. 07757 


-.00607 


0 

10 

.07854 

0. 07854 

0 

0 

.1564 

9 

. 07757 

. 15515 

.00607 

.01214 

.3090 

8 

. 07470 

.14939 

.01154 

.02308 

.4540 

7 

.06998 

. 13996 

. 01589 

. 03177 

.5878 

6 

.06354 

. 12708 

. 01867 

. 03735 

.7071 

5 

. 05554 

.11107 

.01964 

. 03927 

.8090 

4 

.04616 

.09233 

.01867 

. 03735 

.8910 

3 

.03066 

.07131 

. 01589 

.03177 

.9511 

2 

. 02427 

. 04854 

. 01154 

.02308 

.9877 

1 

. 01229 

. 02457 

.00607 

.01214 


= 2 



(13a) 


Wing lift coefficient. — The wing lift coefficient is obtained 
by means of a spanwise integration of the lift distribution, 


where 


7T . mtr 

Vm—^, Sin — 


If the distribution is symmetrical, 


where 




Vms=^V 

The moment of the distribution / can bo found i 

similar manner 

J* ^ J* ^ ^ 


m a 


_l p /2 _A p Cic /2y\ 

2j-^b^{b) 

For asymmetrical lift distributions 

(¥) V. (15a) 

For symmetrical lift distributions 

(¥)»«. (16b) 

m=l \ 0 /m 

Induced-drag coefficient. — The section induced-drag co- 
efficient is equal to the product of the section lift coefficient 
and the induced angle of attack in radians, 


irCiai 

180 

The wing induced-drag coefficient is obtained by means of 
a spanwise integration of the section induced-drag coefficient 
multiplied by the local chord. 



„ 1 TTCiCat 

180 

A p CjC Tra,- 
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For asymmetrical lift distributions 

180 Si { h “V, 

For symmetrical lift distributions 

18ogA6“7. 




(16a) 


(16b) 


Frofilerdrag coefficient. — Tbe section profile-drag coeffi- 
cient can bo obtained from section data for the appropriate 
airfoil section and local Reynolds nmnber. For each span- 
wise station the profile-drag coefficient is read at the section 
hft coefficient previously determined.. The wing profile- 
drag coefficient is then obtained by iheans of a spanwise 
integration of the section profile-drag coefficient multiplied 
by the locah’chord, 



For asymmetrical lift distributions 


For asymmetrical lift distributions 

<7»=S(^) (19a) 

For symmetrical lift distributions 

(19b) 

7»=*1 \ / m 

Rolling-moment coefficient. — The rolling-moment coeffi- 
cient is obtained by means of a spanwise integration 

1 rm 

•4 6 b ^\bj 

= C-r) (20a) 

For an antisymmetrical lift distribution 

C',= -.AS(^) ( 20 b) 

m = l \ O /m 



For symmetrical lift distributions 


(17a) 


(17b) 


Pitching-moment coefficient. — The section pitching-moment 
coefficient about its quarter-chord point can be obtained 
from section data for the appropriate airfoil section and local 
Reynolds number. For each spanwise station the pitching- 
moment coefficient is read at the section lift coefficient 
previously determined and then transferred to the wing 
reference point by the equation 

2/ ' 

, — - [Cj cos (a,— aO+Cdg sin (a,— tti)] ;• 

sin (a«— aj) — cos (a*— ai)] (18) 

where x and z are measured from the wing reference point to 
the 4u^**l®r-chord point of the section under consideration, 
afid upward and backward forces and distances are taken 
as positive. The section pitching-moment coefficient about 
its aerodynamic center may be used instead of in which 
case X and z are measured to the section aerodynamic center. 
The term sin (6:^— a'f) may usually be neglected. The 
wing pitching-moment coefficient is obtained by the spanwise 
integration 



Induced-yawing-moment coefficient. — The induced- 
yawing-moment coefficient is due to the moment of the 
induced-drag distribution, 


r> —A. TTCiCaj 
SbJ_bi2 180 


y dy 


—4 

“4j_i b 180 b \b J 


irA 



( 21 ) 


The induced-yawing-moment coefficient for an antisymmet- 
rical lift distribution is equal to zero and has little meaning 
inasmuch as the lift coefficient is also zero. The induced- 
yawing-moment coefficient is a function of the lift and rolling- 
moment coefficients and must be found for asymmetrical 
lift distributions., 

Prqflle-yawing-moment coefficient. — The profile-yawing- 
moment coefficient is due to the moment of the profile-drag 
distribution, 



APPLICATION OF METHOD USING NONLINEAR SECTION LIFT 
DATA FOR SYMMETRICAL LIFT DISTRIBUTIONS 

The method described is applied herein to a wing, the 
geometric characteristics of which are given in table IV. 
Only symmetrical lift distributions are considered hereinafter 
inasmuch as these are believed to be sufficient for illustrating 


:.ai?256-49— 2 
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TABLE IV.— GEOMETRIC CHARACTERISTICS OF EXAMPLE WING 


Taper ratio, c./ci 

AsiJect ratio. A 

Span, b, ft 

Area, S, sq ft - 

Root chord, c„ ft 

2.5 Root section - 

10. 05 Tip section 

15.00 Geometric twist, ei, deg — ' 

2.143 Edge velocity factor, 

NAOA 4420 

NACA 4412 

-3.50 

- -3.40 

- 1.044 

3. 40X10» 

Mean aerodynamic chord, c', ft 

1. 592 aij , deg 

- - - -3. 90 


2V 

b 

c 

R 

c 

Cj 

c 

6 

c 

c 

£L 

cc' 

ao 

aoc 

'X 

e 

(, deg 

deg 

0 

0.200 

4. 70X10® 

1.0000 

0. 1429 

1.435 

1.932 

0.0969 

0.01385 

0 

0 

0 

.1564 

.195 

4.26 

.9062 

.1295 

1.300 

1.586 

.0973 

.01260 

.0690 

-.24 

-.235 

.3090 

.188 

3. 83 

.8146 

. 1164 

1.169 

1.282 . 

.0978 

.01138 

.1517 

-.53 

-.516 

.4540 

.180 

3. 42 

.7276 

.1040 

1.044 

1. 022 

.0984 

. 01023 

.2496 

-.87 

-.849 

.5878 

. 171 

3.04 

.6473 

.0925 

.929 

.809 

.0991 

.00917 

.3632 

-1. 27 

-1.235 

.7071 

.161 

2.70 

.5757 

.0823 

.826 

.640 

.0999 

.00822 

.4913 

-1. 72 

-1.670 

.8090 

.150 

2. 42 

.5146 

.0736 

.739 

.512 

.1007 

.00740 

.6288 

-2.20 

-2. 138 

.8910 

.139 

2.18 

.4654 

.0665 

.668 

.418 

.1014 

.00674 

. 7658 

-2.68 

-2.604 

.9511 

.129 

2. 02 

.4293 

.0613 

.616 

.356 

.1020 

.00625 

.8862 

-3. 10 

-3. 013 

.9877 

.123 

1.44 

.3061 

.0437 

.439 

.181 

.1021 

.00446 

.9698 

-3.39 

-3. 297 


For tapered wings with straight-line elements from root to construction tip; 


c _/ £i\ 2y / « \ 2y/6 

c. ^ \ C.J b \€,J c. c/c. 

(Alter values of c/c. near tip to allow for rounding.) (Use value of c/c, before rounding tip.) 


the method of calculation. The lift, profile-drag, and pitching- 
moment coefficients for the various wing sections along 
the span were derived from unpublished airfoil data obtained 
in the Langley two-dimensional low-turbulence pressure 
tunnel. The original airfoil data were cross-plotted against 
Reynolds number and thickness ratio inasmuch as both 
varied along the span of the wing. Sample curves are given 
in figures 1 and 2. From these plots the section character- 
istics at the various spanwise stations were determined and 
plotted in the conventional manner. (See fig. 3.) The 
edge-velocity factor E, derived in reference 9 for an elliptic 
wing, has been applied to the section angle of attack for 
each value of section lift coefficient as follows: 

ae=E{oo — a Jq) + a i(, 


LIFT DISTRIBUTION 


Computation of the lift distribution at an angle of attack 
of 3° is shown in table V. This table is designed to be used 
where the multiplication is done by means of a slide rule or 
simple calculating machine. Where calculating machines 
capable of performing accumulative multiplication are avail- 
able, the spaces for the individual products in columns 
to @ may be omitted and the table made smaller. (See 
tables VII and VIII.) The mechanics of computing are 
explained in the table; however, the method for approxi- 
mating the lift-coefficient distribution requires some explana- 
tion. The initially assumed lift-coefficient distribution 
(column d) of first division) can be taken as the distribution 
given by the geometric angles of attack but it is best deter- 
mined by some simple method which will give a close approxi- 
mation to the actual distribution. The initial distribution 
given in table V was approximated by 


^- 1 - 1.8 




where Cj(„) is the lift coefficient read from the section cui’ves 


for the geometric angles of attack. This equation weights 
the lift distribution according to the average of the chord 
distribution of the wing under consideration and that of an 
elliptic wing of the same aspect ratio and span. When 
the lift distributions at several angles of attack are to be 
computed and after they have been obtained for two angles, 
the initially assumed Ci distribution for subsequent angles 
can be more accurately estimated in the following manner: 
Values of down wash angle are first estimated by extra- 
polating from values for the preceding wing angles, and 
then, for the resulting effective angles of attack, the lift 
coefficients are read from the section curves. 

The lift coefficients in column @ of table V, read from 
section lift curves for the effective angles of attack, will 
usually not check the assumed values for the first approxima- 
tion. In order to select assumed values for subsequent ap- 
proximations, the following simple method has been found 
to yield satisfactory results. An incremental value of lift 
coefficient Ac,^ is obtained according to the following relation 

_ (® - ®) m-l + 3<@ - ®) - 0) „,4,1 

where circled numbers represent column numbers in table V 
and where K has the following values at the spanwise 
stations 



K 

0 to 0.8910 

8 to 10 

.9511 

11 to 13 

.9877 

14 to 16 


and (@ — ®)m is tbe difference between the check and as- 
sumed values for the mth spanwise station. The incre- 
mental values so determined are added to the assumed values 
in order to obtain new assumed values to be used in the next 
approximation. This method has been found in practice 
to make the check and assumed values converge in about 
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“> 



(a) Lift. 

(b) Drag. 

(c) Pitching moment. 

Figure l. — Variation of characteristics of NACA 4421 airfoil with Reynolds number. 
(Similar curves plotted for each thickness ratio.) 

three approximations if the first approximation is not too 
much in error. 

WING COEFFICIENTS 

Computations of the wing lift, profile-drag, induced-drag, 
and pitching-moment coeflicients are shown in table VI. 
Since the lateral axis through the wing reference point con- 
tains the quarter-chord points of each section, the x and z 
distances in equation ( 18 ) are zero, and the pitching-moment 



(a) Lift. 

(b) Drag. 

(c) Pitching moment. 

Figure 2. — Variation of characteristics of NACA 44-series airfoil with thickness ratio. 

2y 

R=4.70X10*;-y =0. (Similar curves plotted for Reynolds numbers corresponding to each station.) 

coefficient of the wing is determined solely by the values 
of 

APPLICATION OP METHOD USING LINEAR SECTION LIFT 
DATA FOR SYMMETRICAL LIFT DISTRIBUTIONS 

Although the method described herein was developed 
particularly for use with nonlinear section lift data, it is 


l■lllllnll 
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readily adaptable for use with linear section lift data with a 
resulting reduction in computing time as compared with 
most existing methods. When the section lift curves can be 
assumed linear, it is usually convenient to divide any 
symmetrical lift distribution (as in reference 10) into two 
parts — the additional lift distribution due to angle-of-attack, 
changes and the basic lift distribution due to aerodynamic 
twist. The calculation of these lift distributions is illustrated 


in tables VII to X for the wing, the geometric characteristics 
of which were given in table IV. 

It should be noted that tables VII and VIII are essentially 
the same as table V but are designed primarily for use with 
calculating machines capable of performing accumulative 
multiplication. If such machines are not available, these 
tables may be constructed similar to table V to allow spaces 
for writing the individual products. 



(a) Lift. 

Figure 3.— Section characteristics of example wing. 



(b) • Drag. 

Figure 3.— Cofatinaed. 
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Section pifching-momenf coefficient^ ^'^c/4 


(c) Pitching moment. 
Figure 3.— Concluded. 


LIFT CHARACTERISTICS 


Two lift distributions are required for the determination 
of the additional and basic lift distributions. The first one 
is obtained in table VII for a constant angle of attack 
(e' = 0) and the second one in table VIII for the angle of 
attack distribution due to the aerodynamic twist (aa^=0). 
The check values of cicjh (column @) are obtained by multi- 
plying the effective angle of attack «« by a^clh. The final 


approximations are entered in table IX as ^ and 

It), 


(¥) 


(‘l') 

The 


distribution is the additional lift distribu- 


tion corresponding to a wing lift coeJSicient ^ determined 
in table IX through the use of the multipliers tj^s- It is 


usually convenient to use the additional lift distribution 

Cl c 

— |L. corresponding to a wing lift coefficient of unity. This 
distribution is found by dividing the values of ) 

by <!(...). 

The(%^y , distribution is a combination of the basic 

\ 0 /(‘rO 

lift distribution and an additional lift distribution corre- 
sponding to a wing lift coefficient also determined in 


• Cl c ' 

table IX. The basic lift distribution is then determined 


Of c 

by subtracting the additional lift distribution _ 


from 



Inasmuch as the wing lift curve is assumed to be linear, 
it is defined by its slope and angle of attack for zero lift 
which are also foimd in table IX. The maximum wing lift 
coefficient is estimated according to the method of refer- 
ence 10 which is illustrated in figure 4. The maximum lift 
coefficient is considered to be the wing lift coefficient at 
which some section of the wing becomes the first to reach 
its maximum lift, that is, Cj =Cj . This value of 

Cl is most conveniently determined by finding the minimum 

value of — £!» along the span as illustrated in table IX. 

^'ol 


INDUCED-DRAG COEFFICIENT 


The section induced-drag coefficient is equal to the prod- 
uct of the section lift coefficient and the induced angle of 
attack in radians. The lift distribution for any wing lift 
coefficient is 


CjC 

T 


Cl ,C Ci.C 
■ ‘al /nr I ‘b 

— iT 


(23) 


The corresponding induced angle of attack distribution may 
be written as ; 

^ (24) 


First approximation 


TABLE V— CALCULATION OF LIFT Dll 



. ® 

® 


® 

® 




® 


2 !/ 

b 

a 

Cl 

c 

C|C 



(«»•+«) 

( As - 

sumed ) 

ft 

( TablelV ) 

ft 

(0X0)) 

0 

.1564 

.3090 

.4540 







143.239 

- 58.533 

0 

- 6.950 



3.00 

0.513 

0.1429 

0. 0733 

10.50 

- 4.29 

0 

-.51 







- 1 15.624 

145.025 

- 67.298 

0 


.1564 

2.76 

.517 

.1295 

.0670 

- 7. 75 

9. 72 

- 4.51 

0 








- 64.802 


mmtMl 


.3090 

2.47 

.523 

.1164 

.0609 

0 

- 3.95 

■■ 

- 4.09 







- 12.384 

0 

- 62.917 

160.761 

G 

.4540 

2.13 

.519 

.1040 


-.67 

0 

- 3.40 

8.68 

.O 

"cS 






0 

- 8.320 

0 

- 65.803 

a 

>< 

.5878 

1.73 


.0925 


0 

-.39 

0 

- 3.05 

Ih 

a 





■■I 

- 4.051 

0 

- 7.372 

0 

eO 


1.28 

.477 

.0823 


-.16 

0 

-.29 

0 







0 

- 2.880 

0 

- 7.208 


.8090 

.80 


.0735 


0 

-.09 

0 

-.23 







- 1.638 


- 2.371 

0 


.8910 

.32 

.360 



-.04 

0 


0 





B 


0 

- 1.062 

0 

- 2.016 


.951 1 

-.10 

.281 

IH 


0 


0 

-.03 







-.459 

0 

-.620 

0 


.9877 

-.39 

.228 

.0437 

.0100 


0 

-.01 

0 


1 1.88 

.98 


.TJ 







143.239 

- 58.533 

0 

- 6.950 


0 

3.00 





- 4.17 

0 

-.49 







- 115.624 

145.025 

- 67.298 

BS 


.1564 

2.76 

.516 

.1295 


- 7.72 

9.69 

- 4.50 

BB 







0 

- 64.802 

150.611 

- 67.157 


.3090 

2.47 

.524 

.1164 

.0610 

0 

- 3.95 

9.19 

- 4.10 







- 12.384 

0 

- 62.917 

160.761 

0 

O 

.4540 

2.13 

.517 


.0538 

-.67 

0 

- 3.38 

8.65 

« 

a 


■1 

■ 



0 

- 8.320 

0 

- 65.803 

o 


mgm 



.0463 

0 

-.39 

0 

- 3.05 

a 

a 

Bj 





- 4.051 

0 

- 7.372 

0 

XJ 

a 

p 


1.28 

.478 

.0823 

.0393 

-.16 

0 

-.29 

0 

Q> 

CO 






0 

- 2.880 

0 

- 7.208 


.8090 

.80 

.441 

.0735 

.0324 

0 

-.09 

0 








- 1.638 

0 

- 2.371 

BB 


.8910 

.32 

.382 

■8 


-.04 

0 

-.06 

0 


■i 





0 


^9 



IBD 

-.10 

.292 

.0613 

.0179 

0 

-.02 









-.459 

0 




.9877 

-.39 

.219 

.0437 

.0096 

0 

0 

-.01 

0 

s 

1.61 

1.07 

.95 

.74 


FOR EXAMPLE WING 


bO 


( g ) 


® 




® 

® 

@ 

Xmi 

X ® 





to 

S ®) 


Cheek 

.5878 

.7071 

.8090 

.8910 

.9511 

.9877 

(®-@) 

Cl 

0 

-2.865 

0 

-1.804 

0 

-1.468 




0 

-.21 

0 

-.13 

0 

-.11 

1.88 

1. 12 

0.464 

-10.158 

0 

-4.840 

0 

-3.394 

0 




-.68 

0 

-.32 

0 

-.23 

0 

.98 

1.78 

.531 

0 

-9.916 

0 

-4.968 

0 

-3.768 




0 

-.60 

0 

-.30 

0 

-.23 

.‘90 

1.67 

.522 

-72.472 

0 


0 


BB 




- 3.91 

0 

-.59 

0 


mm 

.77 

1.36 

.514 

177.054 

-82.083 

0 

-13.134 

0 

-7.713 




8.20 

- 3.80 

0 

-.61 

0 

-.36 



.500 

-71.743 

202.571 

-97.965 

0 

-17.388 





- 2.82 

7.96 

- 3.85 

0 

-.68 

■Bll 

.65 

.63 

.474 

0 

-81.434 

243.694 

-125.537 


-26.635 




0 

- 2.57 

7.70 

- 3.97 

0 



.25 

.449 

-7.370 

0 

-96.962 

315.512 

-180.528 





-. 18 

0 

- 2.32 

7.54 

- 4.31 

0 

.42 

-.10 

.413 

0 

-7.599 

0 

-122.880 

463.533 

-329.976 




0 

-.13 

0 

- 2.11 

7.97 

- 5.68 

.77 

-.87 

.326 

-t.49f 

0 

-7.089 

0 

-167.045 

915.651 





0 

-.07 


- 1.67 

9.16 

1.94 

- 2.33 

.165 

■I^B 

.65 

.55 

.42 

.77 

1.94 



0 

2.865 

0 

-1.804 

0 

-1.468 





-.20 

0 

-.13 


-.10 

1.61 

1.39 

.491 

-10.158 

0 

-4.840 

0 

-3.394 

0 




-.68 

0 

-.32 

0 

-.23 

0 

1.07 

1.69 

.523 

0 

-9.916 

0 

-4.968 

0 

WBrn 

mm 




-.60 

0 

-.30 

0 

-.23 


1.52 

.517 

-72.472 

0 

-10.926 

0 

-5.812 

0 

B 



- 3.90 

0 

-.59 

0 

-.31 

0 

H 

1.39 

.517 

177.054 

-82.083 

0 

-13.134 

0 

-7.713 




8.20 

- 3.80 

0 

-.61 

0 

-.36 

.60 

1.13 

.500 

-71.743 

202.571 

-97.965 

0 

-17.388 

BB 

■ 



- 2.82 

7.96 

- 3.85 

0 

-.68 

BS 


.70 

.480 

0 

-81.434 

243.694 

-125.537 

0 

-26.635 

B 



0 

- 2.64 

7.90 

- 4.07 

0 

-.86 

IB 

.19 

.443 

-7.370 

0 

-96.962 

315.512 

-180.528 

0 




-. 19 

0 

- 2.46 

8.01 

- 4.59 

0 

.70 

-.38 

.386 

0 

-7.599 

0 

-122.880 

463.533 

-329.976 




0 

-.14 

0 

- 2.20 

8.30 

- 5.91 

cc 

-.99 

.312 

-1.491 

0 


B|^m 

-167.045 

915.651 



.228 

-.01 

0 

-.07 

0 

- 1.60 

8.79 

1 . 3C 

- 1.72 


58 


61 


70 


1.33 
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Third approximation 



Ci=^S(0X®) =0.490 

C„,-il5M^=.0.0078 


Coo=2(®X®)=0.0077 

C«=S(®X@)=-0.084 


METHOD FOB CALCULATING WING CHARACTERISTICS 




























































































































































































Third approximation Second approximation . First approximation 


TABLE VII— CALCULATION OF LIFT DISTRIBUTION FOR EXAMPLE WING FOR «'=0 



® 

@ 

® 



® 

@ 

® 


® j ® 


® 

® 


® 

® 

® 

21/ 

6 


CiC 

b 

(Assumed) 

Multipliers Xm* 

(S((^X®) 

to 

2(®X@)) 

(‘) 

(®-@) 

Cl 

(ooX®) 

aoc 

b' 

(Table IV) 

Check 

C|C 

b 

(@X@) 

y- 

k y' 

y^ 10 

y' 

9 

8 

7 

6 

5 

4 

3 

2 

1 

y' 

2k / 

Ty 
y 0 
y 

.1564 

.3090 

.4540 

.5878 

.7071 

.8090 

.8910 

.9511 

.9877 


0 

10.000 

0. 1107 

143.239 

-58.533 

0 

-6.950 

0 

-2.865 

0 

-1.804 

0 

-1.468 

2.202 

7.798 

0.7656 

0. 01385 

0. 1080 

.1564 


— 


-115.624 

145.025 

-67.298 

0 

-10.158 

0 

-4.840 


-3.394 

0 

1.475 

8.525 

.8295 

.01260 

.1074 



.0982 

0 

-64.802 

150.61 1 

-67.157 

0 

-9.916 

0 

-4.968 

0 


1.356 

8.644 

.8454 

.01138 

.0984 

.4540 


.0904 

-12.384 

0 

-62.917 

160.761 

-72.472 

0 

-10.926 

0 

-5.812 

BM 

1.236 

8.764 

.8624 

.01023 

.0897 

.5878 



.0819 

0 

-8.320 

0 

-65.803 

177.054 

-82.083 

0 

-13.134 

0 

-7.713 

1.226 

8. 774 

.8695 

.00917 

.0805 

.7071 



.0728 

-4.051 

0 

-7.372 

0 

-71.743 

202.571 

-97.965 

0 

-17.388 

0 

1.257 

8.743 

.8734 

.00822 

.0719 

.8090 



.0632 

0 

-2.880 

0 

-7.208 

0 

-81.434 

243.694 

-125.537 

0 

-26.635 

1.411 

8.589 

.8649 

.00740 

.0636 

.8910 

■ 

1 

.0533 

-1.638 

0 

-2.371 

0 

-7.370 

0 

-96.962 


-180.528 

0’ . 

1.787 

8. 213 

.8328 

.00674 

.0554 

.951 1 

■ 

.0434 

0 

-1.062 

0 

-2.016 

0 

-7.599 

0 


463.533 


3.754 

6.246 

.6371 

.00625 

.0390 

.9877 

B 

■ 


-.459 ! 0 

-.620 ! 0 

-1.491 

0 

-7.089 


-167.045 


.2030 

.00446 

.0089 


0 

■ 

1 


143.239 

-58.533 

0 

-6.950 

0 

-2.865 

0 

-1.804 

0 

-1.468 

2.094 

7.906 

0.7661 

0. 01385 

0.1095 

.1564 



. 1055 

-11 5.624 

145.025 

-67.298 

0 

-10.158 

0 

-4.840 

0 

-3.394 

0 

1. 558 

8.442 

.8214 

.01260 

.1064 

.3090 


— 

.0985 

0 

-64.802 

150.611 

-67.157 

0 

-9.916 

0 

-4.968 

0 

-3.768 

1. 392 

8.608 

.8419 

.01138 

.0980 

.4540 


.0901 

-12.384 

0 

-62.917 

160.761 

-72.472 

0 

-10.926 

0 

-5.812 

0 

1. 213 

8.787 

.8646 

.01023 

.0899 

.5878 


.0813 

0 

-8.320 

0 

-65.803 

177.054 

-82.083 

0 

-13.134 

0 

-7.713 

1.177 

8.823 

.8744 

.00917 

.0809 

.7071 



.0723 

-4.051 

0 

-7.372 

0 

-71.743 

202.571 

-97.965 

0 

-17.388 

0 

1. 205 

8.795 

.8786 

.00822 

.0723 

.8090 



.0634 

0 

-2.880 

0 

-7.208 

0 

-81.434 

243.694 

-125.537 

0 

-26.635 

1. 520 

8.480 

.8539 

.00740 

.0628 

.8910 



.0535 

-1.638 

0 

-2.371 

0 

-7.370 

0 

-96.962 

315.512 

-180.528 

0 

2. 114 

7.886 

.7996 

.00674 

.0532 

.9511 


- 

.0411 

0 

-1.062 

0 

-2.016 

0 

-7.599 

0 

-122.880 

463.533 

-329.976 

3. 379 

6.621 

.6753 

.00625 

.0414 

.9877 


.0232 

-.459 

0 

-.620 

0 

-1.491 

0 

-7.089 

0 

-167.045 

915.651 

4. 832 

5. 168 

.5277 

.00446 

■ .0230 

■ 

0 

E 

■ 


143.239 

-58.533 

0 

-6.950 

0 


0 

-1.804 

0 

-1.468 

2. 060 

7.940 

0. 7694 

0.01385 

0.1100 

.1564 



.1057 

-115.624 

145.025 

-67.298 

0 

-10.158 


-4.840 

0 



1.602 

8.398 

.8171 

.01280 

.1058 

.3090 



.0984 

0 

-64.802 

150.611 

-67.157 

0 

-9.916 

0 

-4.968 


-3.768 

1.377 

8.623 

.8433 

. 01138 

.0981 ' 

.4540 



.0899 

-12.384 



160.761 

-72.472 

0 

-10.926 

0 

-5.812 

0 

1. 203 

8.795 

.8654 

.01023 

.0900 

.5878 




0 


0 

-65.803 

177.054 

-82.083 


HESa 

0 

-7.713 

1.162 

8.838 

.8758 

.00917 

.0810 

.7071 



.0722 



-7.372 

0 

-71.743 

202.571 


■■i 

-17.388 

0 

1. 218 

8.782 

.8773 

. 00822 

.0722 

; .8090 



.0632 


-2.880 



0 

-81.434 


-125.537 

0 

-26.635 

1.492 

8.508 

.8568 

.00740 

.0630 

.8910 



.0534 

-1.638 

0 


■■1 

-7.370 

0 

-96.962 

315.512 

-180.528 

0 

2.111 

7.889 

.7999 

.00674 

.0532 

.9511 


1 .0411 

0 

-1.062 

0 

-2.016 

0 

-7.599 

0 

-122.880 

463.533 

-329.976 

3. 399 

6.601 

.6733 

.00625 

.0413 

• p877 



.0232 

-.459 

0 

-.620 

0 

-1.491 

0 

-7.089 

0 

-167.045 

915.651 

4. 840 

5.160 

.5268 

.00446 

.0230 


|+1.273'»/l-f 

J. j,ClC C \ \ 0 / 

1 list assumed 


!» ,=2®XX„t. 


REPORT NO. 865. ^NATIONAL-: ADVISORY COMMITTEE FOR AERONAUTICS 










































































TABLE Vlir— CALCULATION OF LIFT DISTRIBUTION FOR EXAMPLE WING FOR 



® 

© 

© 

© 

® 

© 

© 

© 

© 




@ 

® 

® 

(1 

® 

® 

6 

t' 

(Table 

IV) 

C|C 

b 

(Assumed) 

Multipliers X™* 

(2((SX®) 

to 

(S®x®)) 

(■) 

ott 

(0-©) 

1 

1 

i 

(aoX 

®) 

OoC 

b 

(Table IV) 

Check 

ClC 

b 

(®X@) 

k 

y 

9 

8 

7 

6 


4 

3 

2 

1 

y 

?y y 

by 

y 0 

.1564 

.3090 

.4540 

.5878 

.7071 

.8090 

.8910 

.9511 

.9877 

j 

First approximation 

0 

0 

0 

143.239 

-58.533 

0 

-6.950 

0 

-2.865 

0 

-1.804 

0 

-1.468 

0.460 

-0.460 

J 

-0.(^46 

0.01385 

-0.0064 

.1564 

—.235 

-.0025 

-115.62 

145.025 

-67.298 

0 

-10.158 

0 

-4.840 

0 

-3.394 

0 

.105 

—.340 

-.0331 

.01200 

-.0043 

.3090 

-.516 

-.0051 

0 

-64.802 

150.61 1 

-67.157 

0 

-9.916 

0 

-4.968 

0 

-3.768 

.012 

-.528 

-.0516 

.01138 

-.0060 

.4540 

-.849 

-.0077 

-12.384 

0 

-62.917 

160.761 

-72.472 

0 

-10.926 

0 

-5.812 

0 

-.107 

-.742 

-.0730 

.01023 

-.0076 

.5878 

-1. 235 

-.0101 

0 

-8.320 

0 

-65.803 

177.054 

-82.083 

0 

-13.134 

0 

-7.713 

-.221 

-I. 014 

-.1005 

.00917 

-.0093 

.7071 

-1. 670 

-.0121 

-4.051 

0 

-7.372 

0 

-71.743 

202.571 

-97.965 

0 

-17.388 

0 

-.373 

-1.297 

-.1296 

.00822 

-.0107 

.8090 

-2.138 

-.0135 

0 

-2.880 

0 

-7.208 

0 

-81.434 

243.694 

-125.537 

0 

-26.635 

-.596 

-1.642 

-.1553 

.00740 

-.0114 

.8910 

-2.604 

-.0139 

-1.638 

0 

-2.371 

0 

-7.370 

0 

-96.962 

315.512 

-180.528 

0 

-.923 

-1.6,81 

-. 1705 

.00674 

-.0113 

.9511 

-3. 013 

-.0131 


-1.062 


-2.016 i 0 1 -7.599 

0 

-122.880 

463.533 

-329.976 

-1. 779 

-1.234 

-.1259 

.00625 

-.0077 

.9877 

-3. 297 

-.0091 

-.459 

0 

-.620 

0 i -1.491 1 0 

-7.089 

0 

-167.045 

915.651 

-3. 653 

.256 

-.0261 

.00146 

.0011 

Second approximation 

0 

0 

-0.0023 

143.239 

-58.533 

0 

-6.950 

0 

-2.865 

0 

-1.804 1 0 

-1.468 

0. 275 

-0.275 

-0.0266 

0.01385 

-0.0038 

.1564 

-.235 

-.0038 

-1 15.624 

145.025 

-67.298 

0 

-10.158 

0 

-4.840 

0 

-3.394 

0 

.075 

-.310 

-.0302 

.01269 

-.0039 

.3090 

-.516 

-.0056 

0 

-64.802 

150.61 1 

-67.157 

0 

-9.916 

0 

-4.968 

0 

-3.768 

.014 

-. 630 

-.0518 

.01138 

-. 0060 

.4540 

-.849 

-.0077 

-12.384 

0 

-62.917 

160.761 

-72.472 

0 

-10.926 

0 

-5.812 

0 

-.095 

-.754 

-.0742 

.01023 

-.0077 

.5878 

-1.235 

-.0097 

, 0 

-8.320 

0 

-65.803 

177.054 

-82.083 

0 

-13.134 

0 

-7.713 

-.202 

-1.033 

-.1024 

.00917 

-.0095 

.7071 

-1.670 

-.0114 

-4.051 

0 

-7.372 

0 

-71.743 

202.571 

-97.965 

0 

-17.388 

0 

-.350 

-1.S20 

-. 1319 

.00822 

-. 0109 

.8090 

-2. 138 

-.0125 

0 

-2.880 

0 

-7.208 

0 

-81.434 

243.694 

-125.537 

0 

-26.635 

-.571 

-1.567 

-. 1578 

.00740 

-.0116 

.8910 

-2.604 

-.0124 

-1.638 

0 

-2.371 

0 

-7.370 

0 

-96.962 

315.512 

-180.528 

0 

-.844 

-1.760 

-.1786 

.00074 

-.0119 

.9511 

-3.013 

-.0109 

0 

- 1 .062 

0 

-2.016 

0 

-7.599 

0 

-122.880 

463.533 

-329.976 

-1.456 

-1.557 

-.1588 

.00525 

-.0097 

.9877 

-3.297 

-.0066 

-.459 


-.620 

0 

-1.491 

0 

-7.089 

0 

-167.045 

915.651 

-2. 014 

-1.283 

-. 1310 

.00446 

-.0057 

Third approximation 

0 

0 

-0.0029 

143.239 

-58.533 

0 

-6.950 

0 

-2.865 

0 

-1.804 

0 

-1.468 

0.210 

-0.210 

-0.0203 

0.01385 

-0.0029 

.1564 

-.235 

-.0040 

-115.624 

145.025 

-67.298 

0 

-10.158 

0 

-4.840 

0 

-3.394 

0 

.085 

-.320 

-.0311 

.01260 

-.0040 

.3090 

-.516 

-.0057 

0 

-64.802 

150.61 1 

-67.157 

0 

-9.916 

0 

-4.968 

0 

-3.768 

.009 

-.525 

-. 0513 

. 01138 

-.0060 

.4540 

-.849 

-.0077 

-12.384 

0 

-62.917 

160.761 

-72.472 

0 

-10.926 

0 

-5.812 

0 

-.095 

-.754 

-.0742 

.01023 

-.0077 

.5878 

-1.235 

-.0096 

0 

-8.320 

0 

-65.803 

177.054 

-82.083 

0 

-13.134 

0 

-7.713 

-.207 

-1.028 

-. 1019 

.00917 

-.0094 

.7071 

-1.670 

-.0111 

-4.051 

0 

-7.372 

0 

-71.743 

202.571 

-97.965 

0 

-17.388 

0 

-.331 

-1.339 

-.1338 

.00822 

-.0110 

.8090 

-2. 138 

-.0121 

0 

-2.880 

0 

-7.208 

0 

-81.434 

243.694 

-125.537 

0 

-26.635 

-.560 

-1.588 

-.1599 

.00740 

-. 0118 

.8910 

-2.604 

-.0120 

-1.638 

0 

-2.371 

0 

-7.370 

0 

-96.962 

315.512 

-180.528 

0 

-.830 

-1.774 

-.1799 

.00674 

-.0120 

.9511 

-3.013 

-.0104 

0 

-1.062 

0 

-2.016 

0 

-7.599 

0 

-122.880 

463.533 

-329.976 

-1.351 

-1.662 

-. 1695 

.00625 

-.0104 

.9877 

-3.297 

-.0063 

-.459 

0 

-.620 

0 

-1.491 

0 

-7.089 

0 

-167.045 

915.651 

-1.915 

-1.382 

-.1411 

.00446 

-.0062 


, -f+1.273Vl-(^y 

First assumed 2A+3.6^ 


*ajj = 2®XXnil;. 
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TABLE IX.— CALCULATION OF LINEAR LIFT CHARACTERISTICS FOR EXAMPLE WING 

[^=10.05; ao, = 10.00;ai^ =-3.90] 


® 



® 

© 


® 

® 

® 

® 

® 

® 

1v 

b 

Multipliers 

T?m. 

\ y 

(Table Vin 

(°‘(-.>) 

(V*),.,, 

(Table VIII) 

(®XCl 

Cl^C 

b 

(®-®) 

c 

1 

(Table IV) 

(1) 

bib 

(1) 

CImax 

(Sectiondata) 

b'ci 

(V) 

0 

0.07854 

0. 1102 

0. 1323 

-0. 0029 

-0. 0105 

0. 0076 

0. 1429 

0.926 

0. 053 

1.421 

1.477 

.1564 

.15515 

. 1057 

. 1269 

-.0040 

-.0100 

.0060 

. 1295 

.980 

.046 

1.418 

1.400 

— 

.3090 

.14939 

.0984 

. 1181 

-.0057 

-. 0093 

. 0036 

.1164 

1.015 

.031 

1.423 

1.371 

.4540 

.13996 

.0899 

. 1079 

-.0077 

-.0085 

.0008 

. 1040 

1.038 

.008 

1.432 

1.372 

.5878 

.12708 

.0811 

.0974 

-.0096 

-.0077 

-.0019 

.0925 

1.053 

—.021 

1.441 

1.388 

.7071 

.11107 

.0722 

.0867 

-.0111 

-.0068 

-. 0043 

.0823 

1.053 

1 

1 1 
1 o 
1 2 
1 
1 

1.436 

1.412 

.8090 

.09233 

. 0632 

.0759 

-.0121 

-. 0060 

-.0081 

.0735 

1.033 

-.083 

1.418 

1.463 

.8910 

.07131 

. 0.634 

.0641 

-.0120 

-.0051 

-.0069 

.0665 

.964 

-. 104 

1.404 

1.564 

.951 1 

.04854 

.0411 

.0493 

-.0104 

-. 0039 

-. 0065 

.0613 

.804 

-.106 

1.419 

1.897 

.9877 

.02457 

.0232 

.0279 

-.0063 

-.0022 

-.0011 

.0437 

.638 

-.094 

1.412 

2.361 

^ =.4S(@X®) =0.833 CL(e,i) =yts(@x®) = -0.079 

C*L («/) 

OCa,, =0.95 

a= — >--=0.0833 •(/..()) a 

F/.„„j=Min. value in @=1.37 “.(i-o) =“'o,+““. [ ;,.oi 


TABLE X.— CALCULATION OF INDUCED-DRAG COEFFICIENT FOR EXAMPLE WING 

[.4 = 10.06; =0.833; =-0.079] 


® 

® 

® 


® 

® 

® 

® 

® 

@ 

® 





“■«i 

“'0,0 

“'■«iCt(.,0 


«'ol« 

ClfcC 

57.3Cd, c 

'al 

57. 3cd . c 

*alb 

57.3cd,. 

'(> 

2y 

Multipliers 


b 

h 

b 

b 

b 

b 


(Table VII) 

(4t;) 

(Table VIII) 

(®xc^o,o) 

(®-®) 

(Table IX) 

(Table IX) 

(®X®) 

((®X®) + 
(®X®)) 

(®X®) 

0 

0.07854 

2. 060 

2. 474 

0. 210 

-0. 195 

0. 405 

0. 1323 

0. 0076 

0. 3273 

0. 0724 

0. 0031 

.1564 

.15515 

1.602 

1.924 

.085 

152 

,237 

. 1269 

.0060 

.2442 

.0416 

.0014 

.3090 

.14939 

1.377 

1.653 

.009 

-. 131 

. 140 

. 1181 

.0036 

. 1952 

.0225 

.0005 

.4540 

.13996 

1.203 

1.445 

-.09.5 

-. 114 

.019 

. 1079 

.0008 

. 1559 

.0032 

0 

.5878 

.12708 

1. 162 

1.395 

-.207 

-.110 

-.097 

. 0974 

-. 0019 

. 1359 

-.0121 

.0002 

.7071 

.11107 

1.218 

1.463 

-.331 

-. 116 

-.215 

.0867 

-. 0042 

.1268 

-.0248 

.0009 

.8090 

.09233 

1.492 

1.792 

-.550 

-.142 

-.408 

.0759 

-.0061 

. 1360 

-.0419 

.0025 

.8910 

.07131 

2. Ill 

2. 535 

-.830 

-.200 

-.630 

.0641 

-.0069 

. 1625 

-.0579 

.0043 

.9511 

.04854 

3.399 

4. 081 

-1.351 

-.322 

-1.029 

.0493 

-. 0065 

. 2012 

-. 0773 

.0067 

.9877 

.02457 

4.840 

5.812 

-1.915 

-.459 

— 1.456 

.0279 

-. 0041 

. 1622 

-.0645 

.0060 


= 0.0322 Ci.2-0.0003 Ct +0.0003 


The values of and 

same manner as and 
distribution is therefore 


, are determined in table X in the 

0 

C\ c 

in table IX. The induced-drag 


c_fi ai 

h ~ h 57.3 


or 


b 



Oz.2 + 


^al b /~y 




(25) 


where 



Cl ,C a 


^al 


b 57.3 


(26) 


and 


b b 57.3~^ b 57.3 


Cd. c 

'■b 


b 57.3 


(27) 

(28) 


The calculation of each of these induced-drag distributions 
is illustrated in table X together with the numerical inte- 
gration of each distribution to obtain the wing induced-drag 
coefficient. 
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PROFILE-DRAG AND PITCHING-MOMENT COEFFICIENTS 

The profile-drag and pitching-moment coefficients for the 
wing depend directly upon the section data and therefore 
their calculation is the same whether linear or nonlinear 
section lift data are used. the linear case the section 

lift coefficient is 

Cl^Ciai^L+Ci^ 

for any wing coefficient (7^. By use of this value for 
the profile-drag and pitching-moment coefficients are found 
as in table VI. 


DISCUSSION 

The characteristics of three wings with symmetrical lift 
distributions have been calculated by use of both nonlinear 
and linear section lift data and are presented in figure 5 
together with experimental results. These data were taken 
from reference 11. The lift curves calculated by use of non- 
linear section lift data are in close agreement with the 
experimental results over the entire range of lift coefficient, 
whereas those calculated by use of linear section lift data are 
in agreement only over the linear portions of the curves as 
would be expected. 

It must be remembered that the methods presented are 
subject to the limitations of lifting-line theory upon which 
the methods are based ; therefore, the close agreement shown 
in figure 5 should not be expected for wings of low aspect 
ratio or large sweep. The use of the edge-velocity factor 
more or less compensates for some of the effects of aspect 
ratio and, ip fact, appears to overcompensate at the larger 
values of aspect ratio as shown in figure 5. 

Additional comparisons of calculated and experimental 
data are given in reference 11 for wings with symmetrical 
lift distributions, but very little comparable data are avail- 
able foi' wings with asymmetrical lift distributions. Such 
data are very desirable in order to determine the reliability 
with which calculated data may be used to predict experi- 
mental wing characteristics. 

Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., December 20, 1946. 



(a) A =8.04; R=4.32X10*; root section, NACA 4416; tip section, NACA 4412. 

Figubf. 5. — Experimental and calculated characteristics of three wings of taper ratio 2..5 and NACA 44-series airfoil sections. 
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Figure 5. — Continued. 


(e) A — 12.06; /?=2.S7X10*; root section, NACA 4424; tip section, NACA 4412. 



Figure 5— Concluded. 
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Pocdtlve direotions of axes and angles (forces and moments) are shown by arrows 


Axis 


Force 
(parallel 
to axis) 
symbol 

Moment about axis 

Angle 

• Velocities 

- 

Designation 

Sym- 

bol 

Designation 

Sym- 

bol 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(comiw- 
nent along 
axis) 

Angular 

Longitudinals 

X 

X 

Rolling 

. 

L 

Y *Z 

Roll 

, 0 
'f' 

tt 

V 

Lateral 

Y- 

Y 

PitehiTigr 

M 

Z *X 

Pitch 

V 

0 

Normal ' 

Z ^ 

Z 

Yawing 

N 

X — 

Yaw 

10 

r 










Absolute coefficients of moment 

/Tf ^ ^ ^ 

(roUmg) (pitching) (yawrbg) 


Angle of set of control surface (relative to neutral 
position), 5. (Indicate surface by proper subscript.) 


4. PROPELLER SYMBOLS 


D 

V 

VlB 

V' 

V, 

T 

Q 


Diameter 
Geometric pitch 
Pitch ratio 
Inflow velocity 
Slipstream velocity 

Thrust, absolute coefficient 
Torque, absolute coefficient Cq— 


Power, absolute coefficient Op= 


C, Speed-power coefficient 






pV^ 

Pn^ 



V 

n 


Efficiency 

Revolutions per second, rps 
Effective helix angle— tan"^|^ 2 ^:^^ 


6. NUMERICAL RELATIONS 


1 hp=76.04 kg-m/s— 660 ft-lb/sec 
1 metric horsepower = 0.9863 hp 
1 mph=0.4470 mps 
1 mps= 2.2369 mph 


1 lb=0.4536 kg 
1^^=2:2046 ib 
1 mi= 1,609.35 m=5,280 ft 
1 m=3.2808 ft 




